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Abstract 

The assumption of a flat prior distribution plays a critical role in the anthropic prediction of 
the cosmological constant. In a previous paper we analytically calculated the distribution for the 
cosmological constant, including the prior and anthropic selection effects, in a large toy "single- 
jump" landscape model. We showed that it is possible for the fractal prior distribution we found 
to behave as an effectively flat distribution in a wide class of landscapes, but only if the single 
jump size is large enough. We extend this work here by investigating a large ( ~ 10^'^'') toy 
"multi-jump" landscape model. The jump sizes range over three orders of magnitude and an 
overall free parameter c determines the absolute size of the jumps. We will show that for "large" 
c the distribution of probabilities of vacua in the anthropic range is effectively flat, and thus the 
successful anthropic prediction is validated. However, we argue that for small c, the distribution 
may not be smooth. 
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I. INTRODUCTION 



A beautiful feature of inflation is that it is generically eternal [l|, 0, Isj giving rise to the 
"multiverse" . Developments in string theory have also led to the complementary world view 
0, M, @] that the fundamental laws of physics admit a vast array of possible solutions. In such 
models there are of order 10^°° different solutions/vacua with various cosmological constants 
0, S 0]- Each vacuum state represents a possible type of bubble universe governed by its 
own low-energy laws of physics. 

This so-called "string theory landscape" of possibilities is expected to have many high- 
energy metastable false vacua which can decay through bubble nucleation(lol. Bub- 
bles of lower-en ergy vacuum can nucleate and expand in the high-energy vacuum background 



and vice versa^ 13|, [ij] . This recycling process will populate the multiverse with bubbles of 
all different types nested one within the other. 

Most of these bubbles will never be home to observers. For example, bubbles with large 
pos itive cosmological constant do not allow for structures such as galaxies or atoms to form 



|15l . |16| . And bubbles with large negative cosmological constant collapse long before life has 
a chance to evolve. However, because of the vastness of the landscape, there will also be 
many bubbles which do provide a suitable environment in which life can flourish. We are 
not surprised to find ourselves in such a fertile bubble. 

In the context of the multiverse, some physical parameters that were once thought of as 
fundamental universal parameters may simply be "local" environmental parameters. The 
most famous example of one such parameter, is that of the cosmological constant. 

The observed value of the cosmological constant A is about 120 orders of magnitude 
smaller than theoretically expected^ 



Ao ~ 10 



-120 



(1) 



A very appealing explanation for this observation assumes that A is an environmental 



parameter which has different values in different parts of the multiverse [15|, [la, ll7|, ll8|, [l 



20, 21, 22 



The probability for a randomly picked observer to measure a given value of A 
can then be expressed as [l7(] 

Pobs(A) (xP(A)nobs(A), (2) 

where -P(A) is the prior distribution or volume fraction of regions with a given value of 
A and nobs(A) is the anthropic/selection factor, which is proportional to the number of 
observers that will evolve per unit volume. If we ignore the variation of other constants, 
then the density of observers is roughly proportional to the fraction of matter clustered in 
large galaxies, nobs (A) /g(A). Using the Press-Schechter approximation |23| for /g(A) we 
can write [l^l 

A, 



nobs (A) ~ erfc 



(3) 



^ However, if the lower-energy vacuum has negative or zero-energy, recycling cannot take place. We will 

call vacua from which new bubbles can nucleate non-terminal, or recyclable vacua, while those which do 

not recycle will be called terminal vacua. 
^ Here and below we use reduced Planck units, Mp = (SttG)^^ = 1, where G is Newton's gravitational 

constant. The theoretical expectation could be as "large" as 10^^^, for example because of supersymmetry. 

However, a discrepancy of 56 orders of magnitude still needs to be explained. 
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where we have normahzed nobs to be 1 for A = and we have parameterized the anthropic 
suppression with a value Ac. For the parameters used in Refs. 2J, |25|, Ac is about 10 times 
the observed value of A, 

Ac ~ 6 X 10-^20_ 

The prior distribution P(A) depends on the unknown details of the fundamental theory 
and on the dynamics of eternal inflation. However, it has been argued 26|, |27| that it should 
be well approximated by a fiat distribution, 



P(A) ^ const. 



(5) 



because the window where nobs(^) is substantially different from zero, is vastly less than the 
expected Planck scale range of variation of A. Any smooth function varying on some large 
characteristic scale will be nearly constant within a relatively tiny interval. Thus from Eq. 

Pobs(A) cx nobs (A). (6) 

Indeed the observed value of A is compatible with the distribution of Eq. ([3]). This 
successful prediction for A depends crucially on the assumption of a fiat volume distribution 
([5]). If, for example, one uses -P(A) oc A instead of ([5]), the 2a prediction would be A/Aq < 
500, giving no satisfactory explanation for why A is so small [28|. 

Given a specific string theory landscape, we would like to be able to calculate the prior 
probability P(A) so that ultimately we can predict the cosmological constant that we should 
expect to observe according to Eq. ([2]). We will assume -P(A) is given by the relative 
bubble abundances of different vacua. Since an eternally inflating multiverse contains an 
infinite number of each type of vacuum allowed in the landscape, it is necessary to use 
some regularization procedure to compute the prior probability distribution. Many such 



regularization procedures, or probability measures, have been proposed [29|, l30|, l3l|, l32|, |33|, 



34, 35 



therein. 



For a more complete and up to date account see Ref. |36| and the references 
Here we will use the pocket-based measure introduced in Refs. 137|, |38i. Refs. 



40 1 computed prior probabilities^ of different vacua in toy models jlHsl and found 
"staggered" distributions ranging over many orders of magnitude. However, to allow for 



40] used models with a relatively small number of vacua and 
See also Refs [ 



42, 43, 44 for related work. 



numerical solution, Refs. [3 
worked onl y in a first-order approximation 

In Ref. 451] we considered a toy model in which bubble abundances were computed as 
though all changes in A were by some fixed amount^ c. In this simple model, we analytically 
studied probability distributions for a realistic number of vacua, N ~ 10^°°. We found that 
when c is around 1, there is a smooth distribution of vacua in the anthropic range, and the 
anthropic prediction of Eq. applies. But when c is smaller by a few orders of magnitude, 
we found that the -P(A) factor, which favors high A vacua, is more important than nobs(A) 
in Eq. ([2]). This implies we should expect to live in a region with large A, and therefore the 
anthropic procedure would not explain the observed small value of A. 



^ Strictly speaking bubble abundances were calculated. 

^ For such a model many degenerate vacua exist with widely separated cosmological constant. We there- 
fore modified the model by artificially perturbing the A of each vacuum, producing a smooth number 
distribution. 
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In this paper we study a more sophisticated model with a range of jump sizes, parame- 
terized by an overall free parameter c. We call it the multi-step model (or MS model). For 
each value of c, the jump sizes range over three orders of magnitude. 

We will show that even for a model with a wide range of jump sizes, the flatness of the 
prior distribution depends on the free parameter c. We find that for "large" c the prior 
distribution is flat but for "small" c there will be some staggering. However it seems as 
though the extent of staggering is reduced in the multi-step model compared to the single- 



step model (or SS model) of Ref. [45| for comparable small values of c and c. 



We also study a simpler "averaged" version of the model, the averaged multi-step model 
(or AMS model), which is more amenable to analytic evaluation and provides some insight 
into the general case. 

The plan of this paper is as follows: We will define our model in section [Til and arrive at 
a general expression for the prior probability distribution in section IIIII We also analytically 
calculate large c limits for the prior probability distribution and present numerical results 
and heuristic arguments for the small c limit in section IIIII 

In section [IV] we will calculate the distribution for the observed A in the large c regime. 
We end the paper with a discussion in section [V| 

In Appendix [A] we will study the averaged multi-step model and calculate the prior 
probabilities in large and small c regimes. In Appendix[B]we will compare expected observed 
probabilities of vacua reached via n and n + 1 jumps for the AMS model. We also include 
an Appendix in which we outline the method used to calculate bubble abundances. 



II. THE MULTI-STEP MODEL 



A. Preliminary outline 

We will study a version of the Arkani-Hamed-Dimopolous-Kachru (ADK) landscape 
model [gI which has J directions, and N = 2^ vacua. We will choose J ~ 1600, so that 
N ~ 10^°°. Each vacuum in this landscape can be specified by a list of numbers {r^i, . . . , 772}, 
where rji = ±1, and the cosmological constant is 

^ = ^ + ^J2v^c^ (7) 

i 

We will take the average cosmological constant A to be in the range (0, c). Each vacuum 
has J neighbors to which it can tunnel by bubble nucleation. Each nucleation event results 
in an increase or decrease of the cosmological constant^ by ci with 1 < i < J = 1600. 



^ For the specific value c « 10^"^, the model has a set of jump sizes ranging roughly between the "small" 
value of 10"'^ and the "large" value of 1. The smaller jump sizes are in the regime where anthropic 
reasoning breaks down for the single step model studied in Ref. 45| , while the larger jump sizes are in 
the regime where anthropic reasoning was found to be valid in Ref. [45|. In the model of Ref. 45| we 
considered 

A = A + i^r;,c (8) 



Thus all jumps had the same size, c. 
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For the model described above vacua exist only with widely separated cosmological con- 
stants, A, A + c. . . , therefore we do not expect any in the anthropic range. Thus we will 
modify the model by artificially perturbing the A of each vacuum to produce a smooth 
number distribution. Vacua originally clustered at A will be spread out over the range from 
to c. This will cover the anthropic range of vacua with A > 0, and so if the density of 
vacua is high enough we will find some anthropic vacua. We will not, however, take account 
of these perturbations in computing probabilities. 

We will only be interested in the vacua near A = 0, which are those that are in the range 
of (0, c) before implementing the perturbation procedure above. One of these will be the 
so-called "dominant" vacuum (defined in Subsection III CI) with A = A*. 



B. Nucleation rates in the multi-jump model 

Vacua with A < are said to be terminal. There are no transitions out of them. Vacua 
with A > are recyclable. If j labels such a vacuum, it may be possible to nucleate bubbles 
of a new vacuum, say i, inside vacuum j. The transition rate Hij for this process is defined 
as the probability per unit time for an observer who is currently in vacuum j to find herself 
in vacuum i. Using the logarithm of the scale factor as our time variable, 

4:71 

K., = T,,—H-\ (9) 

where Tij is the bubble nucleation rate per unit physical spacetime volume (defined in next 
subsection and the same as \ij in [sTj) and 

H, = (A,/3)^/^ (10) 

is the expansion rate in vacuum j. 

Transitions between neighboring vacua, which change one of the integers rja by ±2 can 
occur through bubble nucleation. The bubbles are bounded by thin branes, with tension 
Ta- Transitions with multiple brane nucleation, in which several rji are changed at once, are 
likely to be strongly suppressed [1^1, and we shall disregard them here. 

The bubble nucleation rate Fj, per unit spacetime volume can be expressed as 13 



F,, = A,,exp-^'^ (11) 

with 

Bij = lij - Sj (12) 
Here, lij is the Coleman-DeLuccia instanton action and 

S3 = (13) 

is the background Euclidean action of de Sitter space with expansion rate Hj. 

In the relevant case of a thin-wall bubble, the instanton action lij has been calculated in 
Refs. 0, Il2[ |. It depends on the values of A inside and outside the bubble and on the brane 
tension r. 
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Let us first consider a bubble which changes one rja from rja = +1 to r^^ = —1. The 
resulting change in the cosmological constant is given by 

\AAa\ = ca (14) 

and the exponent in the tunneling rate (fTTj) can be expressed as 

B^, = Bi[^''^-r{x,y). (15) 



j^fiatspace space bouucc actiou, 



277r2 



j-,flatspace 'a 

""^^ = -^jAAJ-r (16) 



The gravitational correction factor r{x,y) is given by [UJ 



2[{1 + xy) - {1 + 2xy + x^)2] 

r[x,y) = J — (17) 

x'^iy'^ ~ 1)(1 + 2x?/ + x^)2 



with the dimensionless parameters 



and 



^ - 7^ (18) 



4|AA,| 



2A 

(19) 



I A A, 



a 



where A is the background value prior to nucleation. 

The brane tension Tj enters the tunneling exponent through the dimensionless parameter 
X (1181) . We will assume that the potentials in our model have the same shape but differ 
by an overall factor so that Vi{(j)) = gfV{(j)). This gives rise to the set of different AAj and 
Tj. However, in this realization of the class of ADK models, the ratio of rf/AAj remains a 
constant for each Vi{(j)). The results we present will correspond to choosing this constant to 
be 1, which is equivalent to having x = 3/4. Although this choice is somewhat ad hoc^, we 
do not expect the main qualitative feature of the results to depend on the precise value of x 

For later convenience, we can rewrite Eq. f[T^ as 



Bai = -—fix,y) (20) 



where we define (see Fig. [T]) 



f(..y).I<i±M^il±^^?^±£!>!l (21) 

(!/-l)(l + 2l!/ + l2)l 



^ According to Ref. 47|, the BPS maximum bound on Ti enforces a: < 1. One can show that small values 
of X lead to enhancement of nucleation rates, which tend to smooth out the distribution. Thus a; = 3/4 
is a reliable choice because it doesn't "wash out" potential problems with staggering. 
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FIG. 1: The factor f as a function of A/AA for the ranges < A/AA < 1 and 1 < A/AA < 100. 

The prefactors Aij in (|TT1) can be estimated'' as 

A, ~ 1. (22) 

If the vacuum {r]i...r]a-i, rja — '^, rja+i----} still has a positive energy density, then an upward 
transition from {r]i...r]a-i,ria - 2,?7„+i....} to {r/i-.. '7a, ^?a+i----} is also possible. The 
corresponding transition rate is characterized by the same instanton action and the same 
prefactor [l3|, and it follows from Eqs. (fTT]) . (fT2|) and (ITOl) that the upward and downward 
nucleation rates are related by 



Fjj — Fjj exp 



-247r^ 



1 1 



(23) 



where Kj > Aj. As expected the transition rate from 77^ = — 1 up to r^^ = +1 is suppressed 
relative to that from r]a = +1 down to rja = —1. The closer we are to Aj = 0, the more 
suppressed are the upward transitions i ^ j relative to the downward ones. 

We will now approximate the dependence of the tunneling exponent Bai on the parameters 
of the model in the following three regimes: 



1. In the limit, A = \AAa\, y = I, 



49 5 

2. For A > |AA„|, we have y > 1, and Eq. ([21]) gives 



9 1 

f{3/A,y = 1) ^ — ^ -. (24) 



I fx Vs [AA 
r{x,y> 1) ^ -= J- ^ —\ — (25) 



V2V y 4 V A 

where we used x = 3/4 in the last step. 
3. For A ^ |AAa|, we have y ~ —1, and Eq. (12T|) gives 



r{x,y^-l)^^^ (26) 



^ Aij is very hard to compute, but using dimensional arguments [49| Aij ~ 0{Hj) ~ 1 for Kj ^ 1. For 
Aj << 1, we expect this result to hold because the tunneling rate remains finite in the limit Aj 0, 
lAAJ ~ 1. 
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For the first two regimes above, r{x,y) < 1, thus the inclusion of gravity decreases the 
tunnehng exponent causing an enhancement of the nucleation rate*. We note that the use 
of the semi-classical approximation is justified when the tunneling action is large enough: 

Bal » 1. 

C. Definition of tlie dominant vacuum 

If we define the total down-tunneling rate for a vacuum j, 

D, = Yl (27) 

Ai<Aj 

then the dominant vacuum, referred to as vacuum *, is defined as that recyclable vacuum 
whose Dj is the smallest. We will call Dj of the dominant vacuum D^. Since bubble 
nucleation rates are suppressed in low-energy vacua, we expect A^, to be fairly small, however 
we would not expect it to be so small as to be in the anthropic range. This can be understood 
as follows: from Eq.'s (I2UI) and fl26l) we see that the transition rates become independent of 
A once we are in the regime A ^ ^A^-. For example, when A ~ lO^^c (still vastly bigger 
than the anthropic range ~ O{10^^^^) for any reasonable value of c) the transition rates 
only depend on AAj, so it is possible for vacuum y with A^; < Aj, < A to have a smaller 
transition rate than vacuum x. 

In Bousso-Polchinski and Arkani-Hamed-Dimopolous-Kachru type landscapes 0, @] it 
can be shown that this vacuum will have no downward transitions to vacua with positive A. 
To see that this is true, imagine that in some direction A^. can jump downward to Aq, > 0. 
Now if we compare to D^, we see that each term contributing to is less than the 
corresponding term (i.e., the transition rate in the same direction) in because Aq, < A^ 
and jump sizes in the same direction are the same. This implies < -D* which contradicts 
our definition of as the vacuum with the smallest sum of downward transition rates. 
Thus for each Kij in D*, A < AAj. 

Recall that when A ^ AA^, the transition rates Kij only depend on AA as seen from 
Eq.'s (!20|) and (!26l) . In this regime, the larger AA the larger Hij. Thus we expect the 
configuration of the dominant vacuum to have rji = —1 for the largest AAj's, so that these 
large transitions are excluded from the sum in Eq. (12 7p . 

Combining this insight with the fact that A^, should be small, we expect the dominant 
vacuum to have a configuration with r]i = {+... H — — } where the first P = 0.71J ~ 1132 
coordinates are -|-'s and the last M = 0.29J ~ 468 coordinates are — 's. These numbers 
are calculated by setting the maximum number of the largest rji to be negative, whilst still 
ensuring that A > 0. 

D. Distribution of vacua 

For a typical realization of the model with J = 1600, the dominant vacuum has + 
coordinates in roughly the first P ~ 1132 directions, and — coordinates in the last M ^ 468 



There are other regimes in which gravity causes suppression of the nucleation rate 10]. However, in this 
paper, all the transitions which will enter in the comparison of probabilities will have A > |AAa| and thus 
are in the graviationally enhanced regime. 
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directions. We also note that the average jump size of the plus and minus coordinates 
respectively are AAp ^ 566c and AAm ~ 1366c. Thus the other vacua with < A < c 
are reached from the dominant vacuum by taking a different number of up jumps and down 
jumps. 

What happens when we take n jumps up from the dominant vacuum? We get a number 
distribution of vacua which is peaked around 

Aji = A. + uAAm (28) 
with a standard deviation which we denote cxn"^''. We find 



cr, 



{up) 



V^135c. (29) 



What happens when we take m jumps down from one of the vacua in the above distri- 
bution? We get another number distribution of vacua which is peaked around 

A^tr^ = Astart - mAAp (30) 

with a standard deviation which we denote cxm""'"'*. We find 

cr^"""") ^ v^327c. (31) 

So overall for a "path" with n up and m down jumps, we get a total distribution which 
is peaked at 

Apelf = A* + nAAM - mAAp (32) 

with 

ajr^ = (33) 

Ultimately, when we come to calculating the probabilities associated with vacua in our 
landscape, we will only be interested in those which are very close to A = (which means 
in our model those that have cosmological constants in the interval < A < c). Thus we 
are especially interested in the special case of 

A'::r - (34) 

which holds when 

m fa ^ '^^Af ^ 2 (■35') 

AAp ^ ^ 

In other words, if you go up n steps, and down m steps, the resultant distribution will 
have its peak around A = if m = an with a = ~ 2.4. 

What happens if , for a given n we go down an — 1 or an + 1 steps? The peaks of 
these distributions will no longer be at A = 0, but if they are broad enough then they may 
still "straddle" the regime of interest < A < c with a significantly high number density 
(see Figl2]). We will consider all vacua within a Icr spread to be statistically significant. It 
turns out that for n ^ 20, the number distribution of vacua which go down m jumps with 
an — 5 ^ m ^ an + 5, straddle the range < A < c within la of their respective peaks. 
This means that there is a statistically meaningful number of vacua in < A < c which are 
reached via all these different "levels" of trajectories. 
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Number of vacua 



Number of vacua 




FIG. 2: On the left we plot six normal distributions with mean (i — l)566c and standard deviation 
(^n°m^^ (see Eq. (j33p ) with n = 20 and m = an — i + 1, and here 1 < i < 6. On the right we plot 
six normal distributions with mean (1 — i)566c and standard deviation an^m^^ with n = 20 and 
m = an + i — 1, and again 1 < i < 6. We take c = 10^^. The heights of these distributions are 
schematic - as we take more downward steps from a starting point, the number distribution grows, 
as indicated by the increasing height of the distributions from right to left in each plot. 

We will classify all vacua in < A < c by two parameters: n and m, where n is the 
minimum number of up jumps required to reach a given vacuum from the dominant vacuum 
and m is the number of down jumps. We will call nm the level of the vacuum. Thus a 
vacuum of level nm differs from the dominant vacuum in n + m coordinates, n of which are 
+ where the dominant vacuum had — , and another m vice versa. 

The total number of vacua of level nm is thus 

Nnm = Nl^P'^ X iV(f°"'") (36) 

where 

^(....f^f^.^^ (37) 

and 



n J n\{M-n)\ 



We approximate these to be smeared over a range (Jnmi where = crn^m''\ so their 
density is 

Pnm ^nml ^nm — ^/ ^nm • ("^^) 

The likelihood that there is no vacuum in a range of size x is exp{—pnmx), thus the 
median A of the lowest-A vacuum is 

Anm = (In 2)/ Pnm = cr„m (In 2) /A^^^ (40) 

We will assume the lowest-A vacuum is at this median position. Above the lowest-A vacuum 
of level nm, there are N^m — 1 more with higher A, with the typical interval in A being Anm- 
For a typical realization it is sufficient to take these vacua as evenly spaced, so that they 
are at 

Km/ = Km + {i-l)Anm (41) 
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TABLE I: h^nm for various values of n and m with c = 0.01. We have highhghted some values to 
draw attention to the fact that paths reached via different levels can have almost the same A^n 
to within a factor of roughly 3 (and similarly for and A^^): In particular when (n,m) = 
(19,2.4n + 4), Anm ~ 3.4 x 10"^^°. Then when {n,m) = {20, an + 1) Anm ~ 9.5 x 10"^^^ and for 
(n,m) = (21, an - 2), Anm ~ 2.8 x lO-^^i^ ^Iso for {n,m) = {22, an - 5), Anm ~ 9 x IQ-'^^^. As 
we went "diagonally" down, Anm decreased by a factor of ~ 3. The same pattern can be seen for 
the "diagonal" immediately above the one just discussed, starting with (n, m) = (19, an + 5) and 



A. 



1.6 X 10-i2i_ 



where I < i < Nnm- 

This model is qualitatively different from the single step model studied in j45|. In the 
single step model, every vacuum in the range (0, c) that was reached via n up and necessarily 
n down jumps (the so-called level n vacua) would have the same number density p„ (and 
also the same prior probability P„). In fact, there is only one way to achieve a density of pn 
and that is to go up and down n jumps. The next level would have pn+i ~ 500p„. 

On the other hand, in the multi-step model, there are multiple classes of trajectories 
which can result in approximately the same densities Pnm- For example for n = 20 and 
m = an+1 ^ 2An + 1, pnm ~ 7.3 x 10^^^ (and Anm ~ 9.5 x lO'^^^) while for n = 21 and 
m = an — 2 we find pnm ~ 2.4 x 10^^° (and Anm ~ 2.8 x 10"^^^). In Table [T]we display where 
we expect the first vacuum Anm to lie for a range of values of n and m. In Fig. [3] we depict 
where we expect some of these Anm vacua to lie, along with second and several subsequent 
vacua for the specific level according to Eq. (HTI) . 

In principle, vacua reached via these different classes will have different probabilities. 
Also, there can even be a spread of probabilities amongst level nm vacua^, depending on 
the paths details. This spread is minimized for large jump sizes. 



III. PROBABILITIES 



We now use the formalism of Refs. 37|, |39j outlined in Appendix O to calculate the 
relative abundances of different vacua in our toy model. 



In contrast to the single-step model which assigns a unique prior probability to every level n vacuum. 
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FIG. 3: Each line represents where we expect to find vacua. The highest Une is at ~ 3.4 x 
10"^^° with n = 19, m = an + A. The next set of highest hnes start at Anm ~ 9.5 x 10~^^^ 
and then are spaced by the corresponding Anm with n = 20, m = an + 1. Also shown are two 
subsets of vacua reached via n = 21, m = an — 2 and n = 19, m = an + 5. Finally, the last and 
most dense set shown is for the first few vacua reached via n = 22, m = an — 5. The height of the 
lines has no meaning here. We simply want to show how different level nm paths result in sets of 
vacua with similar densities. We have set c = 0.01 here. 

The relative abundance of each vacuum (3 is given by a sum over all chains that connect 
it to the dominant vacuum, Eq. (1C2I) . The minimum number of transitions in such a chain 
is roughly n + m. Longer chains can be formed by jumping one way and then later the 
opposite way in the same direction. These chains will have extra suppression factors because 
of the extra jumps, but it remains a possibility that judiciously chosen "extra" jumps might 
enhance the probabilities enough to over compensate for these extra suppression factors. 
We will take the possibility of "extra" jumps into account when we calculate our numerical 
results. For now we shall we forget about them and include only minimum-length chains. 
Furthermore, the paths that maximize the bubble abundances are those that entail first 
making all the up jumps and then following with a sequence of down jumps. The reason is 
that the up-jump suppression factor, Eq. fICip . is least when the starting A for the jump is 
highest. Thus it is best not to jump down until all necessary up-jumps have been made. 

Also, it can be shown that the path that contributes the most to the bubble abundance 
factor of a specific vacuum entails taking up-jumps in order of decreasing jump size while 
the down-jumps should be taken in increasing step sizes. This is because, jumps taken at 
lower A's are more suppressed than those taken at higher A's, so the best strategy is to get 
to high- A's as soon as possible. 

We can reorganize Eq. flC2p . 

E^/3a ^ab ^rs ^tu ^z* f An\ 

Da- D.Dt- D,'" Ds- D.^^'^Dt - D,' " D, - ^ ^ 

The transition rates to the right of the factor Ust in Eq. (1C2p are upward rates, and those 
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to the left are downward rates. Kgt represents the first downward jump after having made n 
upward jumps from the dominant vacuum. 

We will approximate <C Dj, since the transition rates are suppressed for low A vacua. 



Now consider the denominator Dj. 



Dj oc exp 



-247r^ 
AT 



J mm 



H h exp 



-247r2 

"a7 



-r 



J max 



(43) 



where fjmin is calculated from Eq.s fl^ and ffT^ with AA^ in the smallest possible A A 
direction We expect the smallest possible direction A A ~ Ic. 

Similarly Tjmax is calculated with AA^ in the largest possible AA direction. All other 
terms in between have the same form with the rj factors taking on all relevant values as 
AAa increases from the smallest to largest relevant values. 

Using Eq. (ICip . the product of all up-jump suppression factors is 



S 



Ar 



A. 



exp 



1 

A* 



A, 



(44) 



where A^ax is the maximum A reached for a given path. 

The first term in the exponent above and the factor A~^ will contribute to every vacuum 
reached via some path from the dominant vacuum and will thus disappear when we take 
ratios of probabilities. We define 



S 



exp 



247r' 



1 



A, 



(45) 



and call it an up-jump "enhancement factor". 

Consider the last term in Eq. (H2l) which results from the first upward jump 



exp 



-247r2 



1 

A* 



1 

a: 



'dI. 



A. 
A* 



exp 



-247r' 



1 1 

A, A, 
(46) 

Similarly, the product of k/D's for all the upward jumps can be written in terms of 
downward transition rates and an overall suppression factor as 



'D't 



A. 



A, 



exp 



-247r' 



1 

A, 



A. 



(47) 



Now we see that we can call Kij/Dj a branching ratio, since it is the fraction of one 
possible downward transition rate from vacuum Aj divided by the sum of all the possible 
downward transition rates from the same vacuum. 

The down-jumps are similar, except that the first jump down from Amax doesn't have a 
factor of Dt in the denominator (see Eq. fl42p ). giving 



Kir- 



D , - D .. 



15, 



(48) 



Consider vacuum Aj and let us define Jj to be equal to the number of +'s in the configuration describing 
vacuum Aj. There will be Jj different ways to transition down out of Aj to a lower A neighboring vacuum. 
Each one of these possible transitions will have a different factor of f in the relevant transition rate. This 
factor is least when a jump occurs in the smallest AA direction, and we call it fjmin- 
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A. The prior probability of a vacuum of level nm 



The product of branching ratio's for an upward path of n jumps is 



n 

n=l 



Ki£ 



up 



n 

n=l 
n 

n 

n=l 



exp 






(Ef=i exp 


[--.(AA,)]^ 


+ E5=i exp 




) 



■ p 

X^exp 



"aT 



[f,(AA,)-f,(AA^i)] 



^exp 

9=1 



-247r^ 



[r~,(AA,)-f,(AA^,)] 



-1 



(19) 



where P is the number of plus coordinates in the dominant vacuum configuration, and 
we define AA^j to be the change in A for the n'th upward jump. Also, we have introduced 
the notation fj(AAa) = f[x, y) with y = ^ — 1 and x = 3/4 (see Eq's. fll9|2ip ). 

The second sum takes into account that for each upward jump, a coordinate that was 
— becomes a + coordinate, and therefore transitions downward in these directions must be 
included in Dj when calculating the branching ratio. 

Similarly, the product of branching ratio's for a downward path with m jumps is 



m / down 

n(f, 

m=2 ^ J / m 



n 

m=2 
rh—1 



■ P 

X^exp 



1=1 



A ''^ 



[f,(AA,) -f,(AA^)] 



m'=l 
n 

+ ^exp 

9=1 



-247r^ 
-247r2 

"a~ 



[fj(AA„0 -^i(AA^)] 



r,(AA,) -f,-(AA^)] 



n -1 



(50) 



Since the first downward transition does not have a Dj in the denominator (see Eq. ( l48ll ). 
our product of branching ratios starts with the second downward jump, thus m starts at 2 
instead of 1. 

The first sum includes the contribution to Dj from all the original P + coordinates in 
A,^. However, for each jump down on the return path, the next Aj will have one less + 
coordinate and thus the transition rate down in this "lost" direction, must not be included 
in Dj. Thus we subtract the second sum from the first to account for directions which were 
originally +'s but in which we can no longer jump down. 

The third sum takes into account that the downward path jumps take place after n upward 
jumps have changed n coordinates from — to +, and therefore transitions downward in these 
n directions must be included in Dj when calculating the branching ratio. 

Thus the contribution of a path to the prior probability of a vacuum of level nm is 



Pr, 



oc 



n=l 



n(^)"'n 



m=2 



^ \ down 

d) 



with HLi and YIZ=2 )'^°"'" as given by Eq's. and ([50D, and K^t is the first 

downward jump after having reached Amax and is given by 



A. 



A. 



exp 



1 

A, 



1 



A„ 



(51) 
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Kst ^ exp 



-2Att 



2 



A, 



max 



(52) 



Clearly Eq. flSTl) is very complicated. There are Nnm different vacua belonging to level 
nm and each one has in principle a unique prior probability which depends on the details 
of the transitions made, and not just on the number of up and down jumps (as in the single 
step model). We need to stress the fact that for the MS model we have a distribution of 
prior probabilities for level nm vacua, in contrast to the unique prior probability P„ for level 
n vacua in the single-step model. We cannot explicitly calculate all the possible probabilities 
for level nm vacua^^, so we will consider two limiting cases of interest to us. 

It seems reasonable that if c is "large enough", then individual transition rates are high 
and do not differ very much. Thus we expect that the spread in probabilities for level nm 
vacua will be small. If the spread is small enough we can say that level nm vacua have the 
same prior probability of Pnm- Also different orderings to reach the same vacuum should 
become roughly equally weighted. 

On the other hand, if c is "very small", then the suppressed transition rates will differ 
markedly depending on path details and there will be a large spread in the prior probabilities 
of the different vacua within level nm. 

In the next subsection we will try to quantify these heuristic expectations for the large c 
regime. 

The reader may find it useful at this point to turn to Appendix |X] where we study a 
simpler "averaged" version of our multi-step model analytically. 



B. Large c approximation to the prior probability 

For c ~ 10^^ all the branching ratios in Eq's. fH9l) and fl50|) reduce to l/{number of 
possible directions in which downward transitions can be made} (see Fig. 's H] and [5]) and the 
prior becomes^^ 

(l-^A_(AA^=i)) 



^max 



Pnm oc n\m\ Y[{P + n)"' l[{P-{rh-l)+ n)-' A^,, exp 

n=l m=2 

(53) 

where we have dropped both factors in Eq. ( I5T|) which include A* since they will cancel 
when we take ratios of probabilities. Also, we have included the factor of n\m\ to account 
for the different order in which the transitions can be made to get to each level nm vacuum. 



Even for a specific level nm vacuum, there are n!m! different ways to get to this vacuum, and each path 
could be weighted differently. 

The astute reader might notice that in this limit all the transition rates are only slightly suppressed and 
the validity of the semi-classical calculation of the tunneling exponent may be questionable. We wish to 
emphasize that the branching ratios we have used in calculating the prior should still give the correct 
results, even if it is unclear what the exact value of each individual transition rate is. This can be easily 
understood as follows. If we are in a regime of unsuppressed transitions, then all allowed transition rates 
(from a given A) are the same. Thus the branching ratio for a given transition will just be the fraction 
l/{number of possible directions in which downward transitions can be made}, as we have here. 
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FIG. 4: The purpose of these plots is to help us estimate the branching ratios contributing to the 
upward path for "large" c = 10~^. We show three surface plots spanning the possible —Bj = In Kij 
with the following ranges for A: the first plot has 1133 < A/c < 4000, the second has 4000 < 
A/c < 8000 and the third has 8000 < A/c < 30000. All three surfaces have 1133 < AA/c < 1600. 
Notice that along a contour of constant A, there is very little dependence on AA. This shows that 
k/D ~ l/{number of possible directions in which downward transitions can be made}. 




FIG. 5: The purpose of these plots is to help us estimate the branching ratios contributing to 
the downward path for "large" c = 10~^. We show three surface plots spanning the possible 
—Bj = hiKij with the following ranges for A/c and AA/c : the first plot has 450 < A/c < 4000 and 
450 < AA/c < 1132, the second has 4000 < A/c < 8000 and 450 < AA/c < 1132, while the third 
has 8000 < A/c < 30000 and 1 < AA/c < 1132. Again we see that along a contour of constant A, 
the dependence on AA is weak. 



Strictly speaking some paths will be more heavily weighted than others and we should 
replace the factor nlml by some function f{n,m). We assume here that f{n,m) = jSnlml 
where we take /3 to be a constant fraction for any specific value of c. Thus j3 will drop out 
when we evaluate ratios of probabilities for different values of nm. 
Noting that 



n) 



h=l 



and 



Y[{P - {m- 1) +n) 

m=2 

the prior can be written as 



{{P + n)\y 



p\ 



{P + n)\ 

{P + n){P + n-m)\ 
{P + n)\ 



, P\{P + n-m)\{P + n) 
Pnm oc nlml — —2 A„„^exp 



247r2 



(54) 



(55) 



(56) 
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C. Small c regime and the prior probability 

Let's look at Fig's. H] and [5] again. If we were to re-plot these surfaces for smaller c, the 
only thing that would change would be the scale of the height. Since the tunneling exponent 
Bj (X 1/ c, smaller values of c result in larger Bj and thus more suppression. This indicates 
that as c get's smaller there is more of an absolute difference in the amount of suppression 
for jumps occuring in different directions from the same value of A. The net result is that 
different paths will have different prior probabilities even when comparing paths that have 
the same level nm (this spread in probabilities for a given level was not present in the 
single-step and averaged multi-step models). Thus we cannot come up with a closed form 
simplified analytic expression for the prior probability in the small c case. But we will show 
numerically that when c ~ 10~^ vacua of the same level, reached via different paths, start 
to have prior probabilities which differ by several e-foldings. This allows for the possibility 
of a staggered distribution. 

D. Numerical Interlude 

Using numerics we will reinforce our claim that level nm vacua have almost identical 
probabilities even if reached via different subclasses of trajectories for c ~ 10~^. We will 
also reinforce our claim that already for c ~ 10~^ path dependent differences cause a spread 
in probabilities amongst same level nm vacua. We will consider the vacua reached via n 
up jumps and m = an down jumps with n = 19, 20 and 21. We define three subclasses of 
trajectories as follows: 

1. The "evenly spread out" trajectory: Let's consider the n up jumps first. There are 
M = 468 — 've directions to choose from to make our n jumps. Thus we choose a 
trajectory such that within each interval of consecutive M/n coordinates there will be 
one jump. We assume that each jump occurs at the mid-point of the interval. We 
calculate the probabilities by ordering the up jumps from largest to smallest. We do 
the same for the down jumps, only now we order jumps from smallest to largest and 
the size of the interval from which we pick our coordinates for each jump is given by 
P/m, where P is the number of original plus coordinates in the dominant vacuum. 

2. The "extreme" trajectory: Here we jump up in the largest n/2 possible directions, 
and then in the smallest n/2 directions out of the original M coordinates. All the 
up jumps are ordered in decreasing step size. For the downward path we order jumps 
from smallest to largest starting with the m/2 smallest + coordinates, and then ending 
with the m/2 largest + directions. 

3. The "clumped around the average value" trajectory: Again consider the upward path. 
We calculate the average value of AA for the M directions we could go up in. We 
then consider the n/2 coordinates above to the n/2 coordinates below this average, 
and assume that our jumps are made from largest to smallest out of this consecutive 
range of integers. For the downward path we calculate the average AA for the original 
P directions we could go down in. We then consider the m/2 coordinates above to 
the m/2 coordinates below this average, and assume that our jumps were made from 
smallest to largest out of this consecutive range of integers. 
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subclass 


n=19 


n=20 


n=21 


spread 

extreme 

clumped 


454.63 
453.67 
452.57 


475.73 
474.74 
473.68 


496.84 
495.81 
494.78 



TABLE II: For c = 10~^ we find the round trip prior probability and display — In Pnm for the 
three subclasses of trajectories discussed in the text. Results are shown for various values of n 
and m = an. We see that, for a given nm, there is essentially no difference between the different 
trajectories. But the prior is a factor of ~ 10^ times greater for n = 19 than for n = 20, and 
similarly for n = 20 vs n = 21. This is easily understood as follows. The n + 1 trajectory has 
approximately 3 extra jumps over the n trajectory (one for the extra up direction and roughly 2 
extra in the down direction). Each one of these extra jumps contributes a factor of k/D ~ 0(10^^) 
because there are roughly 1000 directions to jump down in. 



subclass 


n=19 


n=20 


n=21 


spread 


545 


566.58 


587 


extreme 


528.9 


549.85 


570 


clumped 


564.2 


585.27 


606 



TABLE III: For c = 10~^ we find the round trip prior probability and display — In Pnm for the 
three subclasses of trajectories discussed in the text. Results are shown for various values of n 
and m = an. We see a large spread in probabilities amongst same level nm vacua which makes it 
impossible for us to use our procedure for comparing observed probabilities in the anthropic range 
as we did in the single-step model, the AMS model discussed in the appendix, and as we will do 
in the large c multi-step case in Section IIV[ 

For each of these trajectories we numerically calculated the prior probability according to 
Eq. (|5T|) for c = 10~^ and c = 10~^ for several values of n. The results are displayed in Tables 
IIT] and nm where we give the — In Pnm for each case. Thus we see that indeed level nm vacua 
have almost identical probabilities even if reached via different subclasses of trajectories for 
c ~ 10~^, whilst for c ~ 10~^ there is already a spread in the prior probabilities amongst 
level nm vacua. 

Let's now take a look at Fig. [H] where we plot —ln{K/D) for each jump in an up and 
down path with {n,m) = (20,48) and c = 10^^. The jumps up were made in order of 
decreasing jump size and the reverse order for the descent back down. The first few jumps 
have the highest values for the "branching ratio" suppression (we distinguish between the 
suppression caused by the branching ratio factors and the overall up-jump suppression which 
depends on how high up in A we have to jump before the downward descent) because they 
are taking place in lower-A backgrounds than the later small jumps in the path The 
small jumps all take place where A is already relatively large, and their suppression effects 



To be absolutely clear, if we made the necessary small jumps first they would lead to much greater 
suppression. But any vacuum is reached via a sum of all possible paths which lead to it, and thus if any 
large jumps are necessary or possible, then it is possible to make the large jumps first and this path to a 
given vacuum will be larger than paths which entail making the small jumps first. 
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FIG. 6: On the left we show — In^n/D) (ie the branching ratio exponential suppression factors) for 
the "spread" up path plotted vs the number of the jump. On the right we show — \n{K/D) for the 
"spread"' down path plotted vs the "reverse" of the number of jumps (for example 10 means the 
tenth to last downward jump and so on). We can see that transitions occurring from lower-A (those 
closer to the y-axis) are more suppressed. Note also that all these value are pretty similar, which 
is to be expected since we are in the large c regime with c = 10~^. The first two up path points are 
at {7.49, 7.30} but are not depicted. Also, not shown are the four points at {7.44, 7.27, 7.19, 7.14} 
closest to the y-axis for the downward path. 
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FIG. 7: Same as previous plots except now c = lO""'. Note that the suppression is 
larger. The first few points are not displayed. For the up path they are {26.71, 13.91}. The 
highest few values for the down path have also been omitted from the display. They are 
{36.40, 17.65, 12.37, 10.21, 9.14, 8.55, ...} for the last jump, second last jump and so on. 

are marginalized. However, we must point out that even if a jump occurs at a high enough 
A such that transition rates are unsuppressed, the branching ratio for the jump is still a 
fraction, k/ D 10~^, so that a path with more steps always pays some penalty. 



IV. DISTRIBUTION FOR THE OBSERVED A IN THE LARGE c REGIME 

We should now be convinced that in the large c regime vacua of level nm have the same 
prior Pnm- We thus have all the parts necessary to calculate the probability of observing 
each value h.nm,i in ^ typical realization of our toy model. These are given by 



19 



The chance that we hve in a world of a given level nm is then given by 

Pohs{nm) OC Pnrn ^ ^obs(Anm,^) (58) 

I 

We will consider two cases^^. 

1. When Knm -C the sum can be approximated by an integral, 

^ . 1 r 3v/^A, ^^NnmK 3v/?A,ln2 
2^n,hs{Kmi) ~ -r— / rfAnobs(A) = — = — = -— (59) 



and we find 

3v^Ac M!P! P!(P + n-m)!(P + n) 2 

(M-n)!(P-m)! ((P + n)-^ '""^ 

■247r2 



X exp 



A, 



(l-fA_(AA^=i)) 



max 



(60) 



where we used Eqs. f l36| [371 [38| [561) . We have not included the terms involving A*, 
which are the same for all Pobs(^'^)- In Eq- fl60l) . Pobs('^^) is a decreasing function of 
n. 



2. When A„m > Ac, Eq. (1581) will be dominated by the first term, and we find 

$^nobs(A„„.,,)^e-(^-/^=)''' (61) 

e 

and 

^247r2 , /A, 



r2!m!P!(P + n — m)!(P + n) 2 
Pobs(?^m) OC — — — —2 A^^^exp 

{{P + n)\) 



(l-fA_(AA^=i)) 



nm 



(62) 

In Eq. (1621) . Pobs(^'^) increases with increasing n while n is small and the last term in 
the exponent is dominant, but it decreases when n is larger and the other terms are 
dominant. 



A. Comparison of probabilities for level n + 1, a{n + 1) + 7 and n, an + 8 + 7 paths 

From Tabledwe see that we are comparing probabilities of vacua which lie adjacent to one 
another along a "bold diagonal" and thus have similar number densities. For example; for 
n = 19, 7 = 1 we would compare probabilities of the two levels (ra, {an+?))+'~^) = (19, al9+4) 
and {n + l,a{n + 1) + 7) = (20, q;20 + 1). We have introduced the bookkeeping integer 
-5 < 7 < 2. 



Recall from the introduction, Ac ^ 6 x 10 parameterizes the anthropic suppression 
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1. When A, 



n,an+3+7 



> Ac and A„+i 



a(n+l)+7 



> Ac, and using Eg. (162 



Pob,(n + l,a(n+l)+7) 



n + 1 + a + 7)! (P + n + 1 - a(n + 1) - 7)! 



Pobsin,an + 3 + -f) (P + n + 1)(P + n) (an + 3 + 7)! (P + n - an - 3 - 7)! 



n+ 1 
X I J exp 



A, 



n,cm+3+7 

X 



2/3 



A 



n+l,o(n+l)+7 



2/3' 



Ac ; 



(63) 



where we have approximated ^ '^^°^('|+^) 
3A„+i_„(„+i)+-y, and we find 



(n±i)2. For n ~ 20, A„, 



an+3+7 



Pob,(n+l,a(n + l)+7) 
Pob^(n, an + 3 + 7) 



n + 1 



(an + a + 7)! (P + n + 1 - a(n + 1) - 7)! 



(P + n + l)(P + n) (an + 3 + 7)! (P 



+ 1 



n 



0.1 exp 



exp 

A, 



A, 



n+l,o(n+l)+7 \ 



2/3' 



Ac 



2/3' 



n 



an 



3-7)! 



(64) 



thus Pobs{n + l,a(r;, + l) + 7) ^ Pobs{n, an + 3 + 7). Therefore, in this regime Pobs 
grows with decreasing A. 

2. When A„+i_„(„+i)+^ < Ac and A„_c,n+3+7 < ^c, we find using Eq. 

Pobsin + l,a{n + 1) + 'y) n + 1 (an + a + 7)! 



Pobs{n,an + 3 + 'j) 



X 



(P + n + 1)(P + n) (an + 3 + 7)! 
n + i y (P + n + l-a(n + l)-7)! A^,c.n+3+7 
n J (P + n - an - 3 - 7)! A„+i,c(„+i)+7 

(65) 



0.3 



Thus we find that several values of nm contribute nearly equally to the total probability. 



The first of these might be dominated by a single A„. 



an+3+75 



but there will be other vacua 



with closer spaced A. These vacua have similar nobs and identical prior probability, so we 
could easily be in any of them. 

Thus when c is large, we recover approximately the original anthropic predictions with 
a smooth prior P(A). There might be an effect due to the discrete nature of the vacua 
associated with the smallest (n, an + 3 + 7) , where Pobs('^, «^ + 3 + 7) has its peak, but this 
effect is small because level (n, an + 3 + 7) does not dominate the probability distribution. 
Instead the probability is divided across many different levels, while only level n, an + 3 + 7 
has the above effect. In other words, we may be more likely to be in the first (n, an + 3 + 7) 
vacuum, but a good fraction of the probability is still spread amongst many other vacua. 
Thus the distribution is effectively flat. 
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B. Comparison of probabilities for level nm and n, (m + 1) paths 



From Table [T] we see that level ra, (m + 1) lie one row above the level nm vacua. These 
vacua have number densities pnm which differ by roughly a factor of 20. We will now evaluate 
the ratio of probabilities for vacua with levels which differ via one extra jump down. 



1. When Knm > Ac and An(m+i) > Ac, and using Eq. flB^ 



Pobs{n{m + 1)) 
Pobsinm) 



m 



1 



m 



exp 



A. 



Ac 



2/3 



A 



?i(m+l) 

aT" 



2/3 



For n ~ 20, m ~ 50, A„m ~ 



P + n 

20A„(„+i), thus Pobs{n{m + 1)) > Pobs{nm 
2. When Knm > Ac and An(m+i) < Ac, using Eq.'s (16011621) we find 

Pobsinim + 1)) m + 1 3v/7rAcln2 



(66) 



Pobsinm) 



exp 



P + n-m 4:An(m+i) 

and therefore Pobsinim + 1)) ^ Pobsinm). 
3. When An^m+i) < Ac and A„m < Ac, using Eq. (160|) . we find 

Pobsinim + 1)) m + 1 A„m 



7 A 




Va 




-^exp 





(67) 



Pobsinm) 



P + n- m An(m+1) 



1. 



(68) 



We can also have vacua reached via m + 2, m + 3, ... jumps down and all these vacua 
contribute nearly equally to the total probabihty. We could live in any one of these. Thus 
since level nm doesn't dominate over level n(m + 1), n(m + 2),... we again see evidence for 
an effectively fiat distribution. 



V. DISCUSSION 

A key ingredient in the anthropic prediction of the cosmological constant is the assump- 
tion of a flat prior distribution. However, the first attempt to calculate this distribution 
for the Bousso-Polchinski and Arkani-Hamed-Dimopolous-Kachru landscape models 3^, 4^ 
revealed a staggered distribution, suggesting a conflict with anthropic predictions. 

These calculations were constrained by computational limitations and revealed only the 
probabilities of a handful of the most probable vacua^^. In Ref. ji^ we went beyond these 
first order perturbative results by studying a simple toy model which allowed analytic cal- 
culation with a large, realistic number of vacua, ~ 10^°°. We found an interesting fractal 
distribution for the prior -P(A). When including anthropic selection effects to determine 
Pofes(A), we found that agreement with observation depends on the only free parameter of 
the model, the jump size c. 

We showed that when c ~ 1, anthropic reasoning does indeed solve the cosmological 
constant problem. Even though the prior distribution has a rich fractal structure, the states 



It is extremely unlikely that any of these vacua should lie in the anthropic range. 
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of interest have similar vacua sufficiently closely spaced to approximate the flat distribution 
well enough to give the usual anthropic results. 

On the other hand, if c is small, of order 10~^, then the agreement with observation was 
found to break down. In this case the volume fraction of universes with large cosmological 
constant is so high that in the overall probability they are greatly preferred even though the 
selection factor disfavors them. 

We emphasize that the primary cause of the massive spread in probabilities in 



which has been referred to as staggering, comes from differences in total up-jump suppression 
factors. This effect was studied analytically in Jumping to a higher Kmax causes any of 
the descendant vacua to be suppressed relative to descendants of lower Kmax paths. Since 
Amax depends on the free parameter c, so does the relative suppression, and thus we can 
always tune c for smooth or staggered distributions. 

In this paper we have studied a multi-step model with a range of jump sizes, |AAj| = ci 
and 1 < ^ < t/, parameterized by an overall free parameter c. (We also studied a simpler 
averaged version of the model, which allowed more progress analytically, especially for the 
small c scenario.) In j45| we conjectured that, for such a model, landscapes with large jumps 
are more likely to give the standard anthropic results, while those with small jumps are likely 
to predict universes unlike ours. We argued that the transition rates between vacua all have 
terms proportional to 1/c in the exponent. Thus when c is small, the transition rates and 
the probabilities of different vacua are very sensitive to the details of the transitions. When 
c is large, all rates are larger and less variable. This conjecture seems to be borne out by 
our study of the multi-step and AMS models here. 

The conclusion for the multi-jump model, that we can find cases where anthropic rea- 
soning is both valid and invalid depending on some free parameter c is analogous to the 
conclusion reached for the single-jump model [45]. However, there are interesting qualita- 
tive differences in the behavior of the prior distributions. Perhaps we could be so bold as to 
turn the anthropic argument around, and use these results to predict what a sensible value 
for c should be? 

We should say a few words about boundary conditions at the regular Planck regime, 
A ~ 0{Mp). In particular, it may be that for the underlying fundamental theory, which 
we assume gives rise to the effective potential we have used in our ADK model, fluxes 
decompactify at the Planck scale and the effective potential asymptotes off to zero. If 
the fluxes do indeed decompactify it would seem like this Planck boundary would act as 
a sink for probability. We have not included these effects in our calculation because all of 
our transitions have taken place within the Planck boundary since even for c = 10~^ and 
n ~ 20, Amax < 320M|, ~ 0.5^-^^ 

One more speculative thought. If you consider a multi-jump model with less of a difference 
between the ratio of the largest to smallest jump in the range (in the current model the 
ratio is 1600 which is pretty big ) this could lead to a smoothing of the prior distribution 
for a reasonably small value of c. If the ratio of maximum and minimum jump sizes were 
constrained to be say within an order of magnitude, it could turn out that we can define a 
unique prior for a given level nm vacuum, in the small c case too. In so doing we would then 
be able to show that the observed distribution smooths out by comparing Pobs for levels nm 
and n, m + 1 as we did in the large c case. 
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APPENDIX A: THE "AVERAGED" MULTI-STEP MODEL 

The dominant vacuum of the multi-step model we have been considering has roughly the 
following structure: 

A,=c{+ + -...-} (Al) 

with P = 0.71J ^ 1132 +'s and M = 0.29J ^ 468 -'s. Jumps down from the plus 
coordinates range in size from Ic < AA < 1132c and jumps up from the minus coordinates 
range in size from 1133c < AA < 1600c. 

The average jump size of the plus and minus coordinates respectively are AAp ^ 566c 
and AAm ~ 1366c. We will define an "averaged" multi-step model (AMS model) which has 
a dominant vacuum with P plus coordinates all having the same jump size AAp and M 
minus coordinates all having the same jump size AAjvf . 

One of the simplifying characteristics of the AMS model is that vacua reached via n up 
jumps and m down jumps from the dominant vacuum all have the same prior probability, 
Pn, which we will now calculate. 



1. The prior probability of the AMS model 

The product of branching ratio's for an upward path of n jumps is 

-247r2 
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where A^ = uAAm + A^. 

The product of branching ratio's for the downward path with m jumps is 



m=2 m=2 



{P — {m — 1)) + n exp 
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and in this expression A^ = uAAm + A^, — (m — l)AAp = A^ax — [m — l)AAp. 
Thus the prior probability of a vacuum of level n is 
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with HLi i^V and U2=2 {^Y""""" as given by Eq's. (jA2D and (jASD, and is the first 
downward jump after having reached Amax and is given by 



Kst ~ exp 



-247r^ 



A, 



■^A_.(AAp) 



(A5) 



Notice that we have included an extra factor of n\m\ to account for all the equally weighted 
paths that lead to the same vacuum because the up jumps and down jumps can be taken in 
any order. 
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2. Large c approximation to the prior probability of the AMS model 

For c ~ 10"^ all the exponential factors in Eq's. flA2|) and flA3|) reduce to 1 and the prior 
becomes 
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where we have dropped the A,, terms appearing in Eq. flA4p . 
Noting that 
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3. Small c approximation to the prior probability of the AMS model 

For c ~ 10~^ the first exponential factor in Eq. flA2l) becomes large and we can ignore 
the second term. In Eq. flASP the second term becomes very small and can be ignored so 
the prior becomes 



P„ oc n\ni\ JJ^ 

n=l 
^max exp 



Pexp 
247r2 



-247r^ 



[f,(AAp)-f,(AA 



M 



-1 m 



fh=2 



(l-fA_(AAp)) 



(AlO) 



where A^ = uAAm + A* above. 
Noting that 
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we finally have 
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APPENDIX B: DISTRIBUTION FOR THE OBSERVED A IN THE AMS MODEL 

We now have all the parts necessary to calculate the probability of observing each value 
A„ f in a typical realization of our toy AMS model. These are given by 

-Pobs(An,£) OC P„nobs(An/) (Bl) 

The chance that we live in a world of a given level n is then given by 

Pobs(^) OC Pn ^ nobs (An/) (B2) 

t 

We will consider two cases. 

1. When A„ <^ A^, the sum can be approximated by an integral, 

J: -obs(A„,) - - dAn^,M = = = (B3) 

where we have used anm = c and therefore N^/c = 1/A„. We dropped the extra 
subscript of m for Nnm and Anm since in the AMS model, once we specify n, m = an 
is specified too. 

2. When A„ > A^., Eq. (lB2p will be dominated by the first term, and we find 

J]nob.(An,.)^e-(^"/^=)''' (B4) 
I 

1. Large c regime 

1. For A„ < Ac, and including Eqs. (136113711381 R9|) . we find 
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We have not included the terms involving A*, which are the same for all Pobs(n)- In 
Eq. flB5l) . Pobs(n) is a decreasing function of n. 

2. When A„ > A^, using Eq. flA9l) . we find 

^247r2 . ..... /A 



;i-fA_(AAp))-(^y 



n!m!P!(P + n — m)!(P + 77.) 2 
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(B6) 

In Eq. (IB6p . Pobsl'^) increases with increasing n while n is small and the last term in 
the exponent is dominant, but it decreases when n is larger and the other terms are 
dominant. 
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The division between regimes occurs when ~ Ac, i.e., 

cln2/A/'„ ~ Ac ~ 6 X 10"^^° 



(B7) 



With c ~ 10~^, we find n ~ 19/20. The dependence on c is weak, with c ~ 10~^ correspond- 
ing to n ~ 19. For n in this range, changing n by one unit changes A„ by a factor of about 
10^. Thus there is at most one n with A„ ~ Ac. 

Comparing Eq.'s flB5l) and Eq. (]B6P for the same n, we see that they differ by a factor of 
AcNn/cexp (A„/Ac)^/'^ ~ 1 if A„ ~ Ac, so there is no big jump due to switching regimes. 

Now let us start with n = 1 and increase n. Certainly with n = 1, A„ 3> Ac by a huge 
factor, we are in the regime of Eq. (lB6p . and Pobs is infinitesimal. As we increase n, Pobs 
increases. Once n is significantly above 1, we can approximate the increase from one step 
to the next as 

Pohsin + 1) f{n + l){a{n + 1))!(P + n + 1 - a{n + 1))! 
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where we have ignored (A„+i/Ac)^'^^ as much less than (A„/Ac)^'''^ and used 
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~ 1.1 for parameters of interest. In the first term f(AAp) << 1 for n ~ 19, so we will ignore 
this contribution. 

Also, since we are in the regime of large c ~ 10^^, for J = 1600, n ~ 19, the first term in 
the exponent is < —1. Thus, for sufficiently small n, the last term in the exponent dominates 
and Pobs('^ + l)/-Pobs(^) ^ 1- There is only an infinitesimal probability that we will be in a 
vacuum of level n, because there are others that are much more probable. As we increase 

Pohs{n) will continue to increase. 

By the time we reach the point where A„ ~ Ac the exponential term becomes insignificant. 
In that case, we switch to the regime of Eq. fIBSp . where Pobs decreases only slowly with 
increasing n. In this regime. 



Pobs(ri + l) _ N„.+i Pn+i / n + r 



Nn Pn 



n 



(BIO) 



which is about 0.41 for parameters of interest. Thus we find that several values of n con- 
tribute nearly equally to the total probability. The first of these might be dominated by a 
single An, but the others will have a large number of closely spaced A. These vacua have 
similar rzobs and identical prior probability, so we could easily be in any of them. 

Thus when c is large, we recover approximately the original anthropic predictions with 
a smooth prior P(A). There might be an effect due to the discrete nature of the vacua 
associated with the smallest n, where Pobs('^) has its peak, but this effect is small because 
level n does not dominate the probability distribution. Instead the probability is divided 
across many different levels, while only level n has the above effect. 
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2. Small c regime 

1. For A„ < Ac, and including Eqs. ([Ml EZl [38l lAT2l) . we find 
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In Eq. ( IBllI) . Pobs(^) is a decreasing function of n. 
2. Wlien A„ > Ac, using Eq. (1A12I) . we find 
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In Eq. ( ]B12I) . Pohs{n) increases with increasing n while n is small and the last term in 
the exponent is dominant, but it decreases when n is larger and the other terms are 
dominant. 

Now let us start with n = 1 and increase n. Again, with n = 1, A„ ^ Ac by a huge 
factor, we are in the regime of Eq. (]B12|) . and Pobs is infinitesimal. As we increase n, Pobs 
increases. Once n is significantly above 1, we can approximate the increase from one step 
to the next as 
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where we have ignored (A„+i/Ac)^^'^ as much less than (A„/Ac)^^^. (For 1 < < 21 the 
prefactor is 0{lQ-^ /IQ-^)). 

For sufficiently small n, the last term in the exponent dominates and Pobs (^^+1) /Pobs (^) 3> 
1. There is only an infinitesimal probability that we will be in a vacuum of level ra, because 
there are others that are much more probable. As we increase n, Pobs('^) will continue to 
increase. What happens next depends on the magnitude of c. 

Let's suppose that 
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for relevant values of n. 

Thus we can find a value of n such that 
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We will now show that for this n, 



Pobs(ri + l) <PobsH, (B16) 

so that we should find ourselves in a vacuum of at most level n. 

It is not clear from Eq. (IBISP whether A„+i/Ac is more or less than 1, so we might need 
to use either Eq. (IBlip or Eq. (]B12p for Pohsin + 1). We will prove the claim using Eq. 
( IBlip . Since this gives a larger value than Eq. (]B12|) . if Eq. (1B16P holds using Eq. (IBlip . it 
will certainly hold using Eq. (lB12p . Thus we will take 
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Now from Eq. (]B14p . we find that 
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Thus unless (A„/Ac)^/^ is extremely close to the upper bound in Eq. (1B15I) . the exponential 
term in Eq. (IBlTp will be infinitesimal, and Eq. (1B16P will follow. If we do have (A„/Ac)^/^ ~ 
247r^/(n^AAjvf), then A„+i/Ac ^ 1. Then the prefactors in Eq. (1B17P are at most about 
10~^ for parameters of interest, and again Eq. flB16P follows. 

Thus we can say with great confidence that we live in a universe with level n or lower. 
From Eq. fIBlSp . we see immediately that we should observe A > Ac, whereas in fact we 



observe Aq ~ O.lAc. If c is significantly smaller than the limit in Eq. flB14p . then we will see 
a very "non-anthropic" universe. 



APPENDIX C: BUBBLE ABUNDANCES 



We review here the procedure for calculating the volume fraction of vacua of a given kind, 
using the "pocket-based measure" formalism of Refs. js^, H . 

Transition rates depend on the details of the landscape. However, if Aj < Aj, the rate of 
the transition upward from i to j is suppressed relative to the inverse, downward transition, 
by a factor which does not depend on the details of the process jl3P^^ 
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Given the entire set of rates Kj, , we can in pri ncip le compute the bubble abundance pa for 



VJIVCXI Vll^ O^U \J± l-CLV^O rt^j, VVt/ l^Oill 111 (^1 l lH^l 1^ 1^ ^^\_/lll|^U.lJ^ Ull^ UU-UIJI^ Cl U U-lH-ldlH^t/ '■^^ 

each vacuum a, following the methods of Refs. |37l . [si. An exact calculation would require 
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We assume only Lee- Weinberg tunnelings and do not consider Farhi-Guth-Guven (FGG) tunnelings 50 1. 
These FGG tunnelings may be faster in upward transition rates, but their interpretation is unclear 51 
and the resulting spacetime cannot be directly handled by the "pocket-based" measure we employ here. 
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diagonalizing an x matrix. But as in Ref. 3^, we can make the approximation that 
all upward transition rates are tiny compared to all downward transition rates from a given 
vacuum. 

Once we have identified the dominant vacuum, the probability for any vacuum a is given 
by (see the next section) 

{Da-D,){D,-D,)---{D,-D,) ^ ^ 

where the sum is taken over all chains of intermediate vacua a,b, . . . ,z that connect the 
vacuum a to the dominant vacuum. 



1. Bubble abundances by perturbation theory 



As shown in Ref. [37[, the calculation of bubble abundances pj reduces to finding the 
smallest eigenvalue q and the corresponding eigenvector s for a huge N x N recycling tran- 
sition matrix R. Bubble abundances are given by 

a a 

where the summation is over all recyclable vacua which can directly tunnel to j. Hg is the 
Hubble expansion rate in vacuum a and we take ^ 1 because q is an exponentially small 
number. 

In a realistic model, we expect to be very large. In the numerical example of Ref. [39 1 
N ~ 10^, while for a realistic string theory landscape we expect ~ 10^°° js], 0,^8, 
0]. Solving for the dominant eigenvector for such huge matrices is numerically impossible. 



However, in Ref. [39|, |53| the eigenvalue problem was solved via perturbation theory, with 
the upward transition rates (see Eq. fICip ) playing the role of small expansion parameters. 
Here we extend this procedure to all orders of perturbation theory. 

We represent our transition matrix as a sum of an unperturbed matrix and a small 
correction, 

R = R(°) + R(^), (C4) 

where R'^^-* contains all the downward transition rates and R*-"*"^ contains all the upward 
transition rates. We will solve for the zero'th order dominant eigensystem {g*^°\s'^°^} from 
R(°^ and then include contributions from R*--*^^ to all orders of perturbation theory. 
If the vacua are arranged in the order of increasing A, so that 

Ai < A2 < . . . < A^, (C5) 

then R(°^ is an upper triangular matrix. Its eigenvalues are simply equal to its diagonal 
elements, 

j<a 

Hence, the magnitude of the smallest zeroth-order eigenvalue is 

=D, = mm{Dg}. (C7) 
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Downward transitions from * will bring us to the negative-A territory of terminal vacua 
Terminal vacua do not belong in the matrix R; hence, Rp^ = for /? 7^ *, and we see 
that the zeroth order eigenvector has a single nonzero component, 




a* • 



(C8) 



Thus, in fact, we could have included only the diagonal elements of R'-"^ and still obtained 
the correct g^^-* and s^^\ 

Now we would like to include the effect of the upward transition rates in the lower trian- 
gular matrix R(°\ to any given order of perturbation theory. To organize the calculation, 
we note that the eigenvalues of an upper triangular matrix are just the diagonal elements, 
and the eigenvectors are given exactly by the perturbation series, which terminates after N 
terms. Thus we will take just the diagonal elements of R'-"^ as our unperturbed matrix, and 
consider the rest of R'-"^ and R*^-*-^ as a perturbation. By summing all terms in the pertur- 
bation series for this perturbation that involve n elements of R*^"*^-*, we find the perturbation 
term of order n in the small upward jump rates. 

The procedure is the same as that involved in finding an eigenstate of the Schrodinger 
equation in rath-order perturbation theory, working in a basis of wavefunctions which diag- 
onalize the unperturbed equation. See for example Eq. (9.1.16) of Ref. Including all 
orders of perturbation theory, we find the result 



where there can be any number of terms in the sum and the vacuum * is not summed over. 
Combining Eqs. ^ and ^ gives Eq. flC2l) . 
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